Cosmological parameter determination from Planck 
and SDSS data in ACHDM cosmologies 



O 
O 

(N 



L.A. Popa^ 

Institute of Space Sciences, Bucharest-Magurele, R-76900, Romania 

C. Burigana and N. Mandolesi 
Istituto TeSRE, Consiglio Nazionale delle Ricerche, Via Gobetti 101, 1-40129 Bologna, Italy 



X5 

D 
tin 

00 i Received ; accepted 

> 
00 

^ I Submitted to ApJ and received, 23 January 2001 

O 

o 

:^ 

o 
U 

> 
X 



^further address: Istituto TeSRE, Consiglio Nazionale delle Ricerche, Via Gobetti 101, 
1-40129 Bologna, Italy 



- 2 - 



ABSTRACT 

We study the complementarity between the cosmological information obtain- 
able with the Planck surveyour and the large scale structure (LSS) redshift 
surveys in ACHDM cosmologies. 

We compute the initial full phase-space neutrino distribution function for 
ACHDM models by using numerical simulations. As initial condition we adopt 
the HDM density fluctuation power spectrum normahzed on the basis of the 
analysis of the local cluster X-ray temperature function and derive the initial 
neutrino phase-space distribution at each spatial wave number k by using the 
Zel'dovich approximation. These initial neutrino phase-space distributions are 
implemented in the CMBFAST code for the integration of the coupled linearized 
Einstein, Boltzmann and fluid equations in fc-space. We find that the relative 
bias between the CMB temperature fluctuations and the underlying matter den- 
sity fluctuation power spectrum in COBE/DMR normalization is given by the 
CDM component normahzed accordingly to the abundance of rich clusters at the 
present time. 

We use the Fisher information matrix approximation to constrain a multi- 
dimensional parametrization of the ACHDM model, by jointly considering CMB 
and large scale structure data according to the Planck and the SDSS experi- 
mental specifications and by taking into account redshift distortions and nonlin- 
ear effects on the matter power spectrum. We found that, although the CMB 
anisotropy and polarization measurements tend to dominate the constraints on 
most of the cosmological parameters, the additional small scale LSS data help to 
break the parameter degeneracies. 

This work has been done in the framework of the Planck LFI activities. 
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Subject headings: Cosmology: cosmic microwave background - large scale 
structure - dark matter - Elementary particles 

1. Introduction 

The simplest hypothesis for the origin of the large-scale structure of the present 
universe is that it is the result of the gravitational instability of small initial density 
perturbations. In this framework, significant improvements on the study of the formation 
and evolution of the cosmological large-scale structure have been achieved in the recent 
years ^. In the same time, new analysis techniques have been developed in order to extract 
more sensitively the long- wavelength portion of the power spectrum (see, e.g., Hamilton 
et al. 1991, Jain, Mo & White 1995, Peacock & Dodds 1996, Ma 1998). The new generation 
of high precision large scale structure (LSS) redshift surveys such as the Sloan Digital Sky 
survey^ (SDSS) and 2dF survey^ are able to measure the total power spectrum of the 
matter density fluctuations with high accuracy for comoving wavelengts A ^ 20 h^^Mpc 
(h=Ho/100 Km s~^Mpc~^ is dimensionless Hubble constant), probing the strong clustering 
regime effects on the power spectrum. 

We are also increasingly able to probe the primordial fluctuations through Cosmic 
Microwave Background (CMB) anisotropy experiments (see, e.g., Tegmark & Zaldarriaga 
2000 for a recent compilation). The Planck surveyor ^ of ESA will observe the microwave 

^http://star- www.dur.ac.uk/virgo/virgo.html 

^http: / / www.sdss.org 

%ttp://mso. anu.edu.au/ colles/2dFCRS 

^http: / / astro.estec.esa.nl/SA-general/Projects/Planck/ 
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sky anisotropies at frequencies between 30 and 900 GHz with FHWM angular resolutions 
between 33' and 5' and sensitivities per FWHM side squared pixel at ~ lO/iK level at least 
for the channels at frequency less than ~ 300 GHz. 

The efforts to measure the CMB anisotropy at sub-degree angular scales and galaxy 
density fluctuations in our local universe can provide independent probes for the structure 
of the universe on similar comoving scales at different cosmological epochs [the sub-degree 
angular scales on the last scattering surface of the CMB correspond to comoving scales of 
~ 50 — 100 h^^ Mpc in the present local universe: an angle 6 degrees subtends a comoving 
distance of 1059{floh)~^Mpc]. The comparison between these two type of measurements 
is an unique way to test the hypothesis of fluctuation growth by gravitational instability; 
when the time dependence of the galaxy linear bias factor (defined as bj = l/crs, where as 
represents the rms mass fluctuations in a sphere of radius R=8 h~^ Mpc) is understood 
the combined data can be also used to constrain the parameters of competing cosmological 
models (see, e.g., Tegmark 1998, Eisenstein, Hu & Tegmark 1998, Eisenstein, Hu & Tegmark 
1999b, Tegmark & Hamilton 2000). 

The combined analysis of the CMB data and the local density fleld fluctuations, as 
suggested by Juszkiewicz, Gorski & Silk (1987), was used to obtain a corresponding CMB 
map as viewed by a distant observer (Bertschinger, Gorsky & Dekel 1990), starting from 
the density fleld as derived by the POTENT procedure (Bertschinger et al. 1990, Dekel 
1994). Also, by using the likehhood analysis of the MARK HI pecuhar velocity data, the 
density power spectrum for a range of parameters within the framework of CDM models 
normalized to COBE/DMR was translated into a range of angular power spectra of the 
CMB anisotropy and compared with the CMB observations (Zaroubi et al. 1996). The 
density fleld reconstructed via Wiener Filter method was then translated into a map of 
AT/T as viewed by a distant observer on his last-scattering surface. 
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It is difficult in the actual experimental context to explain the observations of the 
cosmological large scale structure as well as the CMB anisotropy in the frame of the 
standard Cold Dark Matter (sCDM) model normahzed to COBE/DMR data (Smoot 
et al.. 1992, Wright et al. 1994, Gorski et al. 1994, Bennett et al. 1996). At small scales, 
both the amplitude and the shape of the CMB anisotropy power spectrum predicted by 
this model are inconsistent with the observations of the Large Scale Structure (LSS) of the 
universe as derived by galaxy surveys (e.g., Scott & White 1994; White et al. 1995; Primack 
et al. 1995). The small scale power excess with respect to the large scale power can be 
reduced by the addition of a Hot Dark Matter (HDM) component to the total mass density 
of the universe in form of massive neutrinos (see, e.g., Primack et al. 1995), also motivated 
by the recent experimental results of the atmospheric neutrino oscillation experiments (see, 
e.g., Fukuda et al. 1998, Ambrosio et al. 1998). 

On the other hand, evidences have been accumulated that we live in a low matter 
density universe (see, e.g., Pukugita, Liu & Sugiyama 1999 and the references therein). 
Indications like the Hubble diagram from Type la supernovae (Riess et al. 1998, Perlmutter 
et al. 1998) and the acoustic peak distribution in the CMB anisotropy power spectrum 
(Hancock et al. 1998, Efstathiou et al. 1999) point to a universe dominated by vacuum 
energy, characterized by a cosmological constant A that keeps the universe close to be 
flat. The combined analysis of the CMB anisotropy experiments and Type la supernovae 
observations (Efstathiou et al. 1999) indicates = 0.251q;]^2 = O.GSlois 

(95% confidence errors) respectively for the total matter and the vacuum energy density 
normalized to the critical density, inferring a Hubble constant value Ho=65 Km s^^ Mpc^^. 

Current cosmological constraints on the cosmological parameters obtained by using the 
most recent CMB anisotropy data (Langc ct al. 2000, Balbi et al. 2000) combined with 
Type la supernovae data implies a best fit model close to a fiat ACHDM model having 
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Qrn ~ 0.33, J1a=0.67 and a neutrino density parameter Qi, 0.1, when the priors Ho=65 
Km s~^ Mpc~^ and Q6h^=0.02 are assumed (Tegmark & Zaldarriaga 2000). 

The introduction of HDM component in the form of neutrinos with the mass in the 
eV range supress the growth of fluctuations on all scales below the neutrino free-streaming 
scale k^^ = inGpa'^/ < v >^ (Bond & Sazlay 1980, Hu & Eisenstain 1998, Ma 1999), where 
G is the gravitational constant, a is the cosmological scale factor {a — qq — 1 today), p 
is the density and < v > is the averaged neutrino speed. The magnitude of the power 
suppression is given by AP/P pa —8^l,y/^lrn ~ — 0.087(mi,/eV)(iV;^/r2mh^), where is the 
number of massive neutrino flavours and m,^ is the neutrino mass. The time dependence of 
the free-streaming length implies that neutrinos cluster gravitationally on smaller length 
scales at latter times (see Ma & Bertschinger 1995, Ma 1996, Ma 1999 for discussions on 
the time dependence of the free-streaming scale). In the cosmological models involving 
a HDM contribution to the total energy density, the existence of neutrino free-streaming 
length scale implies that the growth of the density perturbations depends both on time and 
spatial wave number k. 

In the linear perturbation theory, the CMB anisotropy and matter transfer function 
are computed by the integration of coupled and linearized Einstein, Boltzmann and 
fluid equations (Ma & Bertschinger 1995) that describe the time evolution of the metric 
perturbations in the perturbed density field and the time evolution of the density fields 
in the perturbed spacetime for all the relevant particle species (e.g., photons, baryons, 
cold dark matter and massive neutrinos). The unperturbed energy density and pressure of 
massive neutrinos as well as their perturbed energy density and pressure, energy flux and 
shear stress in the /c-space (see Ma & Bertschinger 1995 for their definitions) are specified 
once the full neutrino phase-space distribution function is known. The full massive neutrino 
phase-space distribution depends on the time, on the neutrino positions and momenta and 
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can be represented in the form (Ma & Bertschinger 1995, Ma 1999): 

9s 1 



/(x, q, a) = /o(e)(l + *(x, q, a)), /o(e) 



hi e^l^BT^ + 1 ' 



where q — apis the neutrino comoving momentum, p is the neutrino conjugate momentum, 
e = (o'^ + a^m^)-'^/^ is the neutrino comoving energy, Ty — a~^T^o is the neutrino temperature, 
Tuo = (4/ll)^/^Xyo = 1.927(Tyo/2.7K) K is the present neutrino temperature, T^o being 
the present radiation temperature, Qs is the number of spin for degree of freedom, and 
hp and are the Planck and Boltzmann constants. In the above equation /o(e) is a 
pure Fermi-Dirac distribution that depends only on the neutrino comoving energy and 
^ represents a perturbation term from this distribution depending on time and neutrino 
comoving momentum and position. In the actual version of the CMBFAST code (Seljak 
& Zaldarriaga 1996) the neutrino phase-space distribution, approximated by a pure 
Fermi-Dirac distribution, does not take into account the neutrino position dependence as 
indicated by the above equation. The inclusion of the full neutrino phase-space distribution 
can lead to the perturbations of the quantities related to the computation of the CMB 
angular power spectrum and matter transfer function. 

The full neutrino phase-space distribution function for Cold -|- Hot Dark Matter (CHDM) 
models was computed before (Ma & Bertschinger 1994) by integrating the neutrino geodesic 
equations in the perturbed background spacetime. It was found a positive correlation 
between the neutrino rms velocities and the neutrino density fluctuations at a redshift 
z ~ 15, revealing the contribution of the perturbations to the Fermi-Dirac distribution. 

In this paper we compute the initial full phase-space neutrino distribution function 
for CHDM and ACHDM models by using numerical simulations. We start from the HDM 
density fluctuation power spectrum with the normalization indicated by the analysis of the 
local cluster X-ray temperature function (Eke, Cole & Frenk 1996) and derive the initial 
neutrino phase-space distribution at each spatial wave number k by using the Zel'dovich 
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approximation (Zel'dovich 1970). The neutrino phase-space distributions obtained in 
this way are then implemented in the CMBFAST code as initial neutrino momentum 
distributions for the integration of the coupled linearized Einstein, Boltzmann and fluid 
equations in the k-space and the computation of CMB power spectra and matter transfer 
functions. 

The cosmological parameters of the CHDM and ACHDM cosmological models considered in 
this paper and the initial conditions are presented in Section 2. The numerical simulation 
approach is described in Section 3. The CMB and matter power spectra obtained from 
numerical simulations are discussed and compared with those obtained in the case a 
pure Fermi-Dirac distribution in Section 4. In Section 5 we study the implications of the 
nonlinear bias factor and nonlinear effects of the matter density fluctuation power spectra 
to constrain the main cosmological parameters, when combined CMB and LSS data are 
take into account. Finally, we discuss the results and draw out our conclusions in Section 6. 



2. The initial conditions 

2.1. Cosmological models 

We present the computation of the full phase-space neutrino distribution function for 
two distinct CHDM and ACHDM models specified by the set of cosmological parameters 
given in Table 1, where we report also the parameters of the corresponding CDM (properly 

sCDM and ACDM) models [i.e. the sCDM (ACDM) model with the same parameter of the 
CHDM (ACHDM) model but with ^^'^°^^/^c;dm _ q ^^^^ dark matter density parameter 

^sCDM/ACDM ^ ^CHDM/ACHDM ^ ^CHDM/ACHDMj ^^^^ ^^^^^ ^^^^^^ 

we consider the contribution of one massless and two massive neutrino species while for the 
corresponding CDM models we consider the contribution of three massless neutrino species. 
We assume adiabatic perturbations and the presence of the scalar modes with the spectral 
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index Us — 1, as predicted by the standard inflationary models (Guth & Pi 1981). 

It is usual to describe the shape of the power spectrum of the density perturbations 
with the shape parameter F = Q^/i (Qr/f^r-o)~^^^ (Bardeen et al. 1986), where Qr is 
the energy density parameter for all the relativistic particles and Q^o corresponds to 
the standard model with photons and three massless neutrino species [flro = 1-68130^, 

= 2.3812 X lO-^h-^eiy, 82.7 = K ~ 1.01]. Observations require 0.22 < T < 0.29 

(Peacock & Dodds 1994). The treatment of the cosmological models with HDM component 
requires the introduction of a second shape parameter, T^, = a^/'^fl,^h'^, that characterizes 
the effect of the neutrino free-streaming on the density fluctuation power spectrum (Ma 
1996). Observations at the present time {z ~ 0) require F,^ < 0.021 (Gawiser 2000). The 
values the shape parameters F and F,^ for our cosmological models are also reported in 
Table 1. 

The choice of the cosmological models presented in Table 1 is motivated by observational 
evidencies. The CHDM model with two 2.4 eV neutrinos {rrii, = 92h^fli, eV) was found 
in remarkably agreement with all the available LSS observations if h~ 0.5 (Primack et 
al. 1995). The ACHDM model was found to be the best fit model of the current CMB 
experimental data combined with the Type la supernovae data when priors on Hq and Qf, 
are assumed (Tegmark & Zaldarriaga 2000). 

[] 

2.2. Matter density fluctuation power spectrum and its normalization 

The power spectrum of the matter density fluctuations is defined as (Bardeen et al. 
1996): 

V{k,a)^Ak-mk,ar^^l^y , (1) 
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Fig. 1. — Panels al) and a2): the weighted matter density fluctuation power spectra and 
the HDM component density fluctuation power spectra. The power spectra in the nonhnear 
regime are indicated with dashed hues. All the power spectra are normalized at as and are 
computed at the present time. The experimental data points are derived from the APM 
galaxy survey. The evolution of the weighted matter density fluctuations and of the HDM 
component density fluctuation is shown in panels bl), cl) and b2), c2) respectively. Panels 
bl) and b2): the growth functions at some redshift values as function of the spatial wave 
number k. Panels cl) and c2): the corresponding growth rates at some redshift values as 
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where ^ is a normalization constant, k is the spatial wave number, T{k, a) is the matter 
transfer function and D{k,a) is the growth factor of the density perturbations assuming a 
value DQ{k) = D{k,a = 1) at the present time. We will refer in the following also to the 
growth rate, defined as f{k, a) — dlnD/dlna. 

The dimensionless power spectrum A^(/c,a), defined as the power per logarithmic 
interval in spatial frequency k (the power variance), is given by (Peebles 1980): 

A^(^,a) = ^r»+^P(a,^). (2) 

The total matter density fluctuations contributed by the different matter components 
(baryons, cold dark matter particles, neutrinos) can be obtained by replacing T{k, a) in the 
equation (1) by the the weighted transfer function given by (Ma 1996): 

T{k, a) = ^n{k, a) + ^TcDM{k, a) + ^n(k, a), (3) 

where Th{k,a), TcDM{k,a) and T^{k,a) are the transfer functions of baryons, cold dark 
matter and neutrinos respectively. The growth factor determines the normalization of 
the amphtude of matter density fluctuations relative to the CMB one (Eisenstein, Hu & 
Tegmark 1008). For the sCDM and ACDM models the growth functions depend only 
on time (see Lahav et al. 1991, Carroll, Press & Turner 1992). The presence of neutrino 
free-streaming alters the growth rate and then for CHDM and ACHDM models the growth 
factors and the growth rates depend on the time and the spatial wave number k (Hu & 
Eisenstein 1998). 

It is usual to deflne the normalization of the density fluctuation power spectrum as 
(78(a), the linear rms mass fluctuations in a sphere of radius R=8/i~^Mpc, since the observed 
rms galaxy counts on this scale is about unit, erg (a) is related to the power spectrum 
A^(/c, a) through: 



(78(a) 



,^^AUk,a)W'i.) 





1/2 



(4) 
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where x = kR and W{x) is the window function. We take a top-hat window function of 
radius R=8h~^Mpc in the real space: 

W{x) = 3(sinx - xcosx)/x^. (5) 

It is also possible to set the amplitude of the initial fluctuation spectrum to the measured 
CMB temperature fluctuations on large scales provided by the COBE/DMR experiment 
(Bunn & White 1997). However, the COBE normalization is sensitive to the power 
spectrum at A; ~ 10~^h Mpc~^, that is not relevant for the galaxy clustering. Alternatively 
to the COBE normalization, the abundance of rich clusters can be used to fix the amplitude 
of the initial mass fluctuations close to the quasilinear scale (see also Jenkins et al. 1998 and 
Peacock 2000 for a discussion of the advantages of this kind of normalization). 

We adopt the values of ag obtained by Eke, Cole & Frenk 1996 for CDM models as 
derived from the analysis of the local X-ray temperature function (for as values obtained 
from slightly different analyses and consistent with these values see also White, Efstathiou 
& Frenk 1993, Viana & Liddle 1996, Pen 1997): 

(Tg™ ^ (0.52 ± 0.04)r]-o-52+o.i3Q^ (flat models) (6) 
cr^DM ^ (0.52 ± 0.04)O-°-^^+°-^^- (open models) 

For the cosmological models involving a HDM component, the value of erg is given by (Ma 
1996): 

^CDM ^ ^HDM /r,\ 
(Tg = CTg X (Tg , (7j 

where cr™'^ and a^^^ are the normalizations for the HDM model and the corresponding 
CDM model respectively. For the purpose of this paper we take a^^^ as given by the 
equations (6) and obtained in the COBE normalization, as computed by the 

CMBFAST code. The values of (7g and af^^ for the cosmological models considered here 
are also reported in Table 1. For comparison, we present also the values of (Tg derived when 
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a^^^ is obtained in the COBE normalization. One can see that while for ACHDM models 
the (jg values derived on the basis of these two normalizations are almost identical, they 
differ by a factor of about two in the case of CHDM models. 

We also indicate in Table 1 the values of the free-streaming wave number kfs and of the 
wave number kni at which the nonlinear effects become important (Peacock & Dodds 1996): 

k,, ^ 0.026 (^)'^'QV2hMpc-S (8) 



nl 



((neff + l)/2)! 



l/(neff+3) 



10, 1 ^ dlnPfA;),, , 



R " ' d\nk 

where /cl is the linear wave number and nesiki) is the effective spectral index of the power 
spectrum. 

Figure 1 presents the weighted matter density fluctuation power spectra and the HDM 
component density fluctuation power spectra for the considered cosmological models, with 
the transfer functions given by the CMBFAST code. All the matter power spectra are 
normalized at erg as given by the equations (6)- (7) and are computed at the present time. 

The nonlinear matter density fluctuation power spectra (dashed lines in Figure 1) are 
computed with the formalism of Peacock & Dodds 1996 by using the appropriate linear 
growth rate functions given by the ratio between the power spectrum of each cosmological 
model and the power spectrum of the corresponding CDM model (Ma 1996) calculated with 
the hnear growth factor of the CDM model given by (Peebles 1980): 



9{^m.^K)---^-^ t (9) 



a 2"~io a?H{af' 



where Q.k = 1 — — is the curvature density parameter and H{a) — (O^a"^ + 
Q-kcT^ + ^^a)^''^ is the Hubble expansion rate normahzed to unit at the present time. We 
test the validity of Peacock & Dodds formula for k/h ^ 5 Mpc~^, finding a very good 
agreement with the analytical approximation to the nonlinear power spectra for CHDM 
models obtained by Ma (1998). 
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We present also in Figure 1 the corresponding growth functions and growth rates for the 
considered cosmological models computed at some representative epochs. 

3. N-body simulations 

We obtain the initial neutrino phase-space distributions at each spatial wave number 
k through numerical simulations based on the standard particle-mesh (PM) method 
(Efstathiou & Eastwood 1981, Hockney & Eastwood 1981) usually used to set the initial 
conditions for the nonhnear evolution of the large-scale structure (see Jenkins et al. 1998 
and the references therein). 

The initial neutrino positions and velocities are generated from the HDM matter density 
fluctuation power spectrum with the normalization given by the equation (6)-(7) by using 
Zel'dovich approximation (Zel'dovich 1970). According to the Zel'dovich approximation, 
the perturbed comoving position of a particle r(ro, a) and its peculiar velocity v{fo, a) are 
related to the fluctuations of the density field SpHOMi'f^o, a, k) through: 

f(fo,A;,a) = fo + -D(A;,a)(i(fo) , v{fo,k,a) = D{k,a)d{ro) , (10) 

— * — * ^ 

Vd(fo) = D~ {k,a)5pHDM{ro,k,a), 

where fo is the coordinate corresponding to the unperturbed comoving position and d{fo) 
is the displacement field. In the simulations we use the set of equations (10) to compute 
the perturbed neutrino comoving positions, the neutrino peculiar velocities and the 
displacement of the density fields for each wave number k at the present time. We assign 
to each neutrino a momentum according to the growth function, when the power of each 
mode is randomly selected from a Gaussian distribution with the mean accordingly to the 
power spectrum of the HDM component (Hoffman & Ribak 1991, Ganon & Hoffman 1993, 
Bertschinger 1995), and add a thermal momentum randomly drawn from a Fermi-Dirac 
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distribution: 



/(?o) 




go = Muv/c and M 



(11) 



where v is the neutrino velocity and c is the speed of Ught. We performed simulations 
with 10 X 32^, 10 X 64^ and 10 x 128^ particles. The neutrinos with identical masses are 
randomly placed on 32"^, 64^ and 128^ grids, 10 per grid point, with comoving spacing tq in 
the range (0.5 - 5) h~^Mpc. We verify the convergence of the results with the variation of 
the number of particles and the comoving spacing. The results presented here are obtained 
from ten simulations with 10 x 128^ neutrinos and a comoving spacing of ro=0.5Mpc. This 
large number of particles ensures enough statistics for the computation of the phase-space 
distributions, while this comoving spacing was found to give the minimum variance of the 
likelihood distribution functions presented in Figure 4. In the computation of the set of 
equation (10) we consider only the growing modes, the nonlinear power spectrum up to 
kmax = 6.28 h Mpc~^ and neglect the contribution of the redshift distortions on the power 
spectrum. 

The neutrino momentum field obtained at each wave number k was sampled in fixed 
equispaced points (we use here and normalized to the neutrino total number 

{Npart = 10 X 128^ in the current simulation). 

Figure 2 shows (left panels) the contour plots of the constant particle probabihties 
in the Sq-5pHDM plane, where 5q — qi{k)— < q{k) >, qi{k) being the peculiar neutrino 
momentum and < q{k) > the bulk neutrino momentum of the lattice: 



We also show (right panels) the momentum distribution functions obtained from numerical 
simulations (continuous line) compared with the thermal momentum distribution function 
(dashed line). We find a displacement of the bulk neutrino momentum towards higher 




part 



(12) 
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values for the same wave number when ag increases. Also, for the same value of as, the 
bulk neutrino momentum is displaced towards smaller values with the increasing of the 
wave number, reflecting the /c-dependence of the growth function. The dominant effect is 
given by the variation of the erg value. 

4. The CMB power spectrum from numerical simulations 

The neutrino momentum distributions obtained from numerical simulations at each 
wave number k are used as initial neutrino phase-space distributions in the CMBFAST 
code for the integration of the Boltzmann, Einstein and fluid equations. As usual, the 
neutrino comoving energy is then E — {q'^ + a^Al^)^/^ (Ma & Bertschinger 1995), where q 
is the neutrino comoving momentum. Following the same procedure implemented in the 
CMBFAST code, we compute in synchronous gauge the perturbations of the energy density, 
pressure, energy flux and shear stress, truncating the Boltzmann hierarchies for massive 
neutrinos at l^ax — 50 for every value of q, obtaining a relative accuracy in the estimation 
of the CMB power spectrum better than ~ 10'^ (Popa et al. 2000). 

Figure 3 presents the evolution with the scale factor of the energy density parameters 
of different components and of their density perturbations, and of the neutrino shear stress 
and energy flux obtained for the CHDM model for the neutrino phase-space distribution 
function modelled as a Fermi-Dirac distribution (continuous lines) or obtained from the 
numerical simulations when the normalization of the matter power spectrum is erg = 0.3 
(dashed lines) or as = 0.81 (dotted-dashed lines). One can see that while the time evolution 
of the various energy density components obtained from numerical simulations are not 
changed, their perturbations as well as the neutrino shear stress and energy flux differ from 
those obtained in the case of a pure Fermi-Dirac distribution. As it was mentioned before, 
the full neutrino phase-space distribution depends also on the neutrino position, leading to 
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k=0.02 MpC' ff,=0.3 




Fig. 2. — Left panels: the contour plots of constant particle probabilities in the 5q — SpnoM 
plane (see also the text). From exterior to interior the contours correspond to: 0.75, 0.5 
and 0.25 probability. Right panels: the momentum distribution functions obtained from 
numerical simulations (continuous line) and the thermal momentum distribution function 
(dashed line). These plots refer to the case of the CHDM cosmological model. 
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Fig. 3. — Panel a): the time evolution of the energy density of the different components (see 
also the text). Panel b): the time evolution of energy density perturbations for different 
components when the neutrino phase-space distribution is a Fermi-Dirac distribution (con- 
tinuous line) or it is obtained from numerical simulations with cxg = 0.3 (dashed line) or with 
(78 = 0.81 (dotted-dashed line). Panel c): the time evolution of the neutrino shear stress 
when the neutrino phase-space distribution is a Fermi-Dirac distribution (continuous line) 
or it is obtained from numerical simulations with erg = 0.3 (dashed line) or with erg = 0.81 
(dotted-dashed line). Panel d): the same as in the panel c) but for neutrino energy flux. 
The distributions presented in panels b) c) and d) are obtained at k — kfg. All the plots 
refer to the CHDM cosmological model. 
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perturbations of the pure Fermi-Dirac distribution, f{x, q, t) — /o(5)[l + ^(^, q, t], where 
fo{q) is the pure Fermi-Dirac distribution. This perturbations are reflected by the time 
evolution of the quantities presented in Figure 3. 

Figure 4 presents the CMB anisotropy power spectra obtained from numerical 
simulations for different normalizations of the matter density fluctuations power spectra, 
for the CHDM (panel al)) and the ACHDM models (panel bl)). Each power spectrum 
is obtained by averaging the power spectra from ten simulations with 10 x 128^ particles 
and normalized to COBE/DMR (Bunn & White 1997). We also show in Figure 4 the 
likelihood dependence on erg for the CHDM model (panel a2)) and ACHDM model (panel 
b2)). For each case the target power spectrum is the corresponding power spectrum given 
by the CMBFAST code normalized to the COBE/DMR data. We found at 1 - o" level the 
following values for as'. 

as = 0.34 ± 0.09 for CHDM 
(78 = 0.3 ±0.07 for ACHDM. 

These values of erg show that the CMB power spectrum in COBE normahzation is well 
recovered when the matter density fluctuation power spectrum is normalized to the cluster 
abundancy data, as indicated in the Table 1. The results presented in Figure 4 show that 
the relative amplitude between the CMB temperature fluctuations and the matter density 
fluctuations normalized to COBE/DMR is well deflned by the CDM component with the 
normalization indicated by the abundance of the rich clusters at the present time. This 
result confirms the predictions of the structure formation theories in which the CDM driven 
by the adiabatic fluctuations leads to the formation of the CMB anisotropy and large scale 
structure (see, e.g., Dodelson, Gates & Turner 1996). Figure 5 presents the time evolution 
of the bias factor bj = l/crg (panel a)) and of the parameter Q^f /6/ (panel b)), usually 
measured from the peculiar velocity data, obtained for the CHDM and the ACHDM model. 
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Fig. 4. — The CMB anisotropy power spectra obtained from numerical simulations for 
CHDM model (panel al)) and ACHDM model (panel a2)) for different values of ag. From 
the top to the bottom (at the first Doppler peak): as— 0.2, 0.3, 0.5. The Ukelihood depen- 
dence on as obtained in CHDM model (panel bl)) and ACHDM model (panel b2)). All the 
reported power spectra are obtained averaging over ten simulations and are normalized to 
COBE/DMR data. 
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For both models the normahzation at the present time is given by the equations (6)-(7). 

5. Cosmological parameter determination from CMB and LSS data in 

ACHDM cosmologies 

The results obtained in the previous section show that the relative bias between the 
CMB anisotropy power spectrum and the matter density fluctuation power spectrum in 
the COBE/DMR normalization is given by the CDM component with the normalization 
obtained from the abundance of the reach clusters at the present time. This confirms the 
idea to jointly use the cosmological information contained in the CMB power spectrum and 
in the matter power spectrum to constrain degenerated sets of cosmological parameters. 

There are few factors that alter the relative amplitude between CMB temperature 
fiuctuations and the underlying mass density fiuctuations: the growth rate of perturbations 
depends upon and and implicitly upon Q^- (see equation (9)). The degeneracy 
between Q\ and Qk leeds to an indeterminacy of the Hubble constant (Efstathiou & Bond 
1999): 

h= {Um + '^A + i^ky^^, (13) 

where cuj = Qih"^ are the physical densities of the relevant components. The CMB anisotropy 
measurements can constrain — 00^ + ujc + and cUb from the morphology of the Doppler 
peaks. 

On the other hand, the observed power spectrum in the redshift space differs from the 
theoretical power spectrum of mass fiuctuations because of the nonlinear evolution, reshift- 
space mapping and the bias. For the purpose of this paper we consider a scale-independent 
bias factor 6/ = l/ug, and compute the effects of the nonlinear evolution by using the 
Peacock & Dodds formula (see Section 2). Following Feldman, Kaiser & Peacock (1994) we 
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Tablc 1: Parameters of the considered cosmological models and other relevant quantities 



Model Qb h T cr™^ as kfs Ki 

(Mpc-i) (Mpc-i) 



CHDM 0.05 0.75 0.2 0.5 0.58 0.05 0.68 0.36±0.02 0.81 0.06 0.46 

ACHDM 0.03 0.2 0.1 0.67 0.65 0.25 0.04 0.33 0.32±0.02 0.31 0.03 0.35 

sCDM 0.05 0.95 0.5 0.5 - - 0.52±0.03 1.17 - 0.46 

ACDM 0.03 0.3 0. 0.67 0.65 0.21 - - 0.88±0.06 0.92 - 0.35 




Fig. 5. — The redshift dependence of the bias factor bj (panel a)) and of the parameter Vt^^ /hi 
(panel b)) obtained for the CHDM and the ACHDM models, when the CDM component at 
the present time is normalized to the cluster abundancy data (see also the text). 
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correct the nonlinear power spectrum of mass fluctuations for redshift distortions. On large 
scales the redshift distortions cause a linear increase of the power given by: 

On small scales the redshift distortions cause: 

^W-Wf^, (16) 

where a measures the spatial rms distortion and erfc is the complementary error function. 
The typical values of a, indicated by the rms small scale velocity measurements (Feldman, 
Kaiser & Peacock 1994), is of few h~^Mpc. The indeterminacy of h causes the movement of 
the matter density power spectrum in the redshift space, changing the normalization scale, 
the neutrino free-streaming scale k^, and the scale k^i at which the nonlinear effects become 
important (see Section 2). 

The relative normalization between the CMB and matter power spectra depends also on 
cUm- While cUm and cuj, can be constrained from the CMB Doppler peak distribution, the 

strong dependence of the matter power spectrum on at intermediate and small scales 
constrains uJc- 

The spectral index of the scalar modes Ug, the spectral index of the tensorial modes rit, 
the ratio T/S between tensorial and scalar contributions as well as the reionization effects 
(an optical depth to the last scattering r 7^ ) affect also the relative bias between the 
CMB temperature fluctuations and the matter density fluctuations. Within the CMB 
polarization data rit, T/S and r can well be constrained, while the indeterminacy of most 
of the cosmological parameters is reduced (see, e.g., Zaldarriaga, Spergel & Seljak 1997, 
Zaldarriaga & Seljak 1997). 

We compute the precision at which the fundamental cosmological parameters could 
be extracted from the combined CMB anisotropy and LSS data given by Planck and 



36/ bh] 



(14) 
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SDSS, taking into account the dependence of the CMB and matter power spectrum on 
all the factors enumerated before. The previous papers on the parameter estimation from 
CMB and LSS data (see, e.g., Scott et al. 1995, Bond, Jaffe & Knox 1998, Lineweaver 1998, 
Webster et al. 1998, Wang et al. 1998, Eisenstein, Hu & Tegmark 1999b) in general consider 
the normahzation to the COBE data, neglecting the cosmological contribution to the 
estimation of the cosmological parameters of the nonlinear effects and redshift distortions. 

We use the Fisher matrix approximation to compute the errors on the estimates of 
the cosmological parameters (see, e.g., Efstathiou & Bond 1999, Popa et al. 1999) for the 
experimental specifications of Planck and SDSS. In this case, the Fisher matrix elements 
are the sum of two terms, accounting for the Planck measurements of the CMB anisotropy 
and polarization and the galaxy power spectrum as derived from SDSS (Tegmark 1997): 



The minimum error that can be obtained on a parameter Sj when we need to determine all 
parameters jointly, is given by: 



depending on the experimental specifications, the target model and the class of considered 
cosmological models. If both anisotropy and polarization power spectra arc used, the first 
Fisher information matrix term is given by (Zaldarriaga & Seljak 1997, Zaldarriaga 1997): 



where X and Y stands for T, C and B power spectra and Cov~^ is the inverse of 
the covariance matrix. The Planck experimental parameters and the procedure used to 
compute the Fisher matrix term given by the equation (18) arc described in Popa et al. 
(2000); in particular, we consider here, for simplicity, the Planck "cosmological channels' 
only (between 70 and 217 GHz) and neglect the foreground contamination. 



-1 _ zTiPlanck , pSDSS 



(16) 




(17) 




(18) 
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The second term of the Fisher information matrix in equation (16) is given by (Tegmark 
1997): 

r^SDSS o /"^""•^ 91nP(A;) ,,,ainP(fc)^, , 

i7;S.DSS ^ 27r / ^ ^ w{k) ^ M ln/c, (19) 

.Ik (J Si USj 

where kmin, kmax are the minimum and maximum wave number used to compute the matter 
density fluctuation power spectra (we take a fixed kmin = lO^'^hMpc^^ while the value of 
kmax can vary) and w{k) is the selection function for the Bright Red Galaxy (BRG) sample 
of the SDSS. According to Tegmark (1997): 

= (20) 

Here Veff{k) is the effective volume of the BRG sample used for measuring the power at 
the wave number k corresponding to the wavelength A = 27i/k: 

where n{r) is the selection function of the survey that gives the expectation value of 
the number density of galaxies. We consider the BRG sample of SDSS volume-limited 

at lOOOh^^Mpc, containing 10^ galaxies (Tegmark 1997) with a bias factor given by the 
equations (6)-(7) and n(r) is the expectation value of the Poisson distribution. 

The fiducial model is the flat ACHDM model with: fi^ = 0.047, 9.^ = 0.233, = 0.67, 
fli, = 0.05 {nil, — 2.1eV), Hq = 65 Km s~^ Mpc~^, one relativistic neutrino flavor, Nrei = 1, 
and two massive neutrino flavors, F = 0.25, Fj, = 0.021, as — 0.42, kfg ~ O.OSMpc"^ and 
kni ~ 0.32Mpc'~^. We assume primordial adiabatic perturbations, the presence of the 
scalar modes with the scalar spectral index Ug = 1, the contribution of the tensorial modes 
(gravitational waves) with the spectral index rit — —0.09, a tensor-to-scalar ratio T/S=0.1 
and an optical depth to the last scattering t — 0. The independent parameters in our 
computation are: uib, Uc, uiu, f-^A, ng, rit, T/S, r and a. We fix ag at the value given by 
the equations (6)- (7) for each corresponding cosmological model and consider a linear bias 

bi = l/cTg. 



n{rP{k) 
1 + n(r)V(k) 



dV, 



(21) 
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Fig. 6. — Panel a): the CMB anisotropy power spectrum of the fiducial model (continuous 
line) and its decomposition into scalar and tensorial contributions (dashed lines), compared 
with the CMB anisotropy experimental data from COBE, Boomerang and MAXIMAL Panel 
b): the polarization power spectra of the fiducial model for E-channel (continuous line) 
and cross-correlation C-channel (dashed and dotted-dashed lines; the dottcd-dashcd line 
represents the negative part of the power spectrum). Panel c): the matter fluctuations 
power spectrum of the fiducial model and Veff{k)/X^ obtained for the fiducial model and the 
experimental specifications of the BRG sample of SDSS. The sohd curves are for the hnear 
regime, the dashed curves are obtained for the nonlinear regime and the dotted-dashed 
curves are obtained when both nonlinear effects and redshift distortions are considered, for 
few choices of u; from top to bottom (7=0,0.2,0.5,1 h~^Mpc (see also the text). The power 
spectra are computed at the present time and normalized according to the equations (6)-(7). 
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Panel a) of Figure 6 presents the CMB anisotropy power spectrum of the fiducial model 
(continuous line) and its decomposition into scalar and tensorial contributions (the dashed 
lines) compared with the CMB anisotropy experimental data from COBE, Boomerang and 
MAXIMAL Panel b) shows the polarization power spectra of the fiducial cosmological 
model. Panel c) of Figure 6 presents the matter fluctuations power spectrum of the flducial 
model and the effective volume of the BRG sample of SDSS obtained for the fiducial model. 
The solid curves are for the linear regime, the dashed curves are obtained for the nonlinear 
regime and the dash-dotted curves are obtained when both nonlinear effects and redshift 
distortions are considered. For the last case we plot the matter power spectrum and Ve///A^ 
for few choices of parameter a in the equation (15). We take for the fiducial model cr ~ 0. 
This makes the correction factor in equation (15) of order unit for the fiducial model. 

Figure 7 presents the 1 — a errors on the estimates of the cosmological parameters 
that can be obtained by SDSS alone as a function of kmax- The sohd curves are obtained 
assuming the linear matter power spectra in equations (19)-(21). The dashed curves are 
obtained when the nonlinear effects and the redshift distortions on the matter power spectra 
are taken into account. 

For the computation of the Fisher matrix elements given by the equation (19) we take 
two-sided derivatives of V{k) considering a step size of Asj a; 5%. 

Parameters as nt, T/S and r, not considered in Figure 7, do not affect the matter power 
spectrum for the SDSS experimental specifications considered here. The results presented in 
Figure 7 show that the SDSS power spectrum information alone can give modest constraints 
on most of the cosmological parameters in both linear and nonlinear regime, strongly 
dependent on the choice of k^ax value. 
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Fig. 7. — 1 — a errors on the estimates of the cosmological parameters that can be obtain by 
SDSS alone from the hnear matter power spectrum (continuous hnes) and from the matter 
power spectrum corrected for nonhnear effects and redshift distortions (dashed hnes) as 
function of k^ax- The maximum k accessible to SDSS (~ 0.3 Mpc~^) is indicated by the 
arrow. 
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Fig. 8. — 1 — a errors on the estimates of the cosmological parameters from the combined 
data sets of Planck and SDSS when: the hnear matter power spectrum is used (continuous 
hnes), the matter power spectrum is corrected for nonhnear effects and redshift distortions 
(dashed hnes). In each case we indicate by the filled circle the 1 — a error obtained by 
Planck alone from anisotropy and polarization data. The arrow indicates the maximum k 
accessible to the SDSS. 
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Figure 8 presents the 1 — a errors on the estimates of the cosmological parameters that 
can be obtained by Planck alone when the CMB anisotiopy and polarization data are 
taken into account and when the combined data sets (CMB anisotropy and polarization and 
matter density fluctuations) are taken into account for different values of kmax- We consider 
for the Planck surveyor a sky coverage fsky — 0.85. Figure 8 shows that, although the 
CMB anisotropy and polarization measurements tend to dominate the constraints on most 
of the cosmological parameters, the additional small scale LSS data {k > 0.1 Mpc~^) help 
to break the parameter degeneracies. 

6. Conclusions 

In this paper we study the complementarity between the cosmological information 

obtainable with the Planck surveyour and the large scale structure redshift surveys in 
ACHDM cosmologies. 

We compute the initial full phase-space neutrino distribution function for CHDM 
and ACHDM models by using numerical simulations. We start from the HDM density 
fluctuation power spectrum in the normalization indicated by the analysis of the local cluster 
X-ray temperature function and derive the initial neutrino phase-space distribution at 
each spatial wave number k by using Zel'dovich approximation. The neutrino phase-space 
distributions obtained in this way are implemented in the CMBFAST code as initial 
neutrino momentum distributions for the integration of the coupled linearized Einstein, 
Boltzmann and fluid equations in /c-space and the computation of CMB power spectra and 
matter transfer functions. 

Wc flnd that the relative bias between the CMB anisotropy power spectrum and the 
matter density fluctuations power spectrum in the COBE/DMR normalization is given by 
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the CDM component normalized to the abundancy of the rich clusters at the present time. 

Taking into account the redshift distortions and nonlinear evolutionary effects 
on the matter density fluctuations power spectrum, we constrain an 11-dimensional 
parametrization of the ACHDM model, when the combined CMB and LSS data are taken 
into account in the Planck and BRG sample of SDSS experimental specifications. We find 
that, depending on the maximum spatial wave number, the combined CMB and LSS data 
can better constrain most of the cosmological parameters. 

It is a pleasure to thank U. Seljak and M. Zaldarriaga for the use of the CMBFAST 
code (v3.2) employed in the computation of the CMB power spectra and of the matter 
transfer functions. 
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